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The upcoming p + Pb run at the LHC will probe the nuclear gluon distribution at very small 
1 Bjorken x (from x ~ 1CP 4 at mid-rapidity down to x ~ 1CP 6 in the proton fragmentation region) 

, and will allow to test approaches based on parton saturation. Here, we present the predictions 

of the KLN model for hadron multiplicities and multiplicity distributions in p + Pb collisions at a 
center-of-mass energy of 4.4 TeV. We also compare the model to the existing pp, dA and AA data 
O I from RHIC and LHC. 

Ov 

Very soon, the Large Hadron Collider will record the first data on p + Pb collisions at the center-of-mass energy of 
4.4 TeV. This data will allow to probe the nuclear gluon distributions at very small Bjorken x: from x ~ 10~ 4 at mid- 
' (— | ' rapidity down to x ~ 10~ 6 in the proton fragmentation region. Since the QCD evolution makes parton distributions 
O j increase at small x, the LHC experiments will allow to probe the nuclear wave functions at unprecedented parton 
' . densities. These measurements are crucial for testing the current theoretical approaches to high energy QCD. 
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Due to the breaking of scale invariance by quantum effects, QCD possesses a dimensionful scale Aqcd that deter- 
mines the characteristic distance ~ Aqq D at which the dynamics becomes non-perturbative. The asymptotic freedom 
[HQ makes the perturbative expansion valid only if a hard external scale Q 2 Aq CD is present. Multiparticle pro- 
duction in hadron collisions is dominated by soft interactions and so in general is not amenable to the weak coupling 
treatment. However when the density of partons in the transverse plane Q 2 becomes large compared to Aq CD , it 
£f~) . regularizes the infrared behavior of the parton transverse momentum distributions at the "saturation momentum" Q 
Q s and thus prevents the running coupling of QCD from growing large, a s (Q s ) = g 2 /4:ir <C 1 0-[5|. The gluon field 
C") . A in this weak coupling regime has a large occupation number, A ~ 1/g > 1 and can be treated as a classical "Color 
Glass Condensate" (CGC) [10]. 
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While the complete theory of multi-particle production based on the ideas outlined above is still being developed, 
its main ingredients are clear and can serve as the basis for phenomenology. This was the motivation for the KLN 
model fa-flG] combining the Glauber approach to proton-nucleus and nucleus-nucleus collisions (for a complete set of 
formulae see e.g. (Tlj ) with a simple ansatz for the unintegrated parton distributions that accounts for the existence 
of a new dimensionful scale - the saturation momentum. The KLN model was successful in describing the RHIC 
data fL2Hl5| on the centrality and rapidity dependence of charged hadron production in heavy ion collisions. The 
predictions for Pb Pb and p Pb collisions at the LHC were made in [l(| . The comparison to the first LHC data [I?} on 
hadron production in Pb-Pb collisions revealed that while the KLN model describes the centrality dependence rather 
well, the overall normalization exceeds the observed one by about 10-15 %. This implies that the energy dependence 
of the saturation momentum assumed in fl6j was slightly too steep 1 . 

Regarding pA collisions, we also have to remember that the number of "participants" (the nucleons that underwent 
at least one inelastic interaction) in this case is much smaller than in A A collisions, and that fluctuations are much 
more important. Therefore a Monte-Carlo (MC) based formulation [l8[ of the numerical integration of the KLN 
model [lj| can be expected to provide more accurate predictions. Indeed, the MC method leads to a better agreement 
between the data and the model prediction [2(| in d Au collisions at RHIC. The MC based KLN model [l8| has been 
used to generate initial conditions for the hydrodynamical description of collective flow, see e.g. (21H23j . 



1 While it is evident that the model has to be refined, let us put this discrepancy in perspective by noting that some of the early pre- RHIC 
predictions for the LHC that did not take into account the concepts of parton saturation and coherence overestimated the measured 
hadron multiplicity by almost an order of magnitude. 
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The goal of this letter is to provide updated predictions for p Pb collisions at the LHC. Let us explicitly list the 
differences between the present and the previous [l6| papers: i) we consider the c.m.s. energy of the forthcoming 
p Pb run - 4.4 TeV; ii) we have reduced the intercept describing the energy dependence of saturation momentum 
by ~ 20%; iii) we employ the MC method of evaluating the number of participants. Of course, after making these 
changes we have to make sure that the RHIC data is still adequately described - therefore we present the comparison 
to the RHIC AA and dA as well. While these changes may seem insignificant, the p Pb LHC data present a chance 
to test saturation ideas, and this requires quantitative predictions made to the best of our current knowledge. 

Let us briefly recall the basic ingredients of the KLN approach; for details, see 

HG1. The multiplicity per unit 

rapidity 

dN K f^2fr, da \ K 47riY c a s f°° dp 2 



dy~S J dPt { E jr p )=J -Nf^l J -pf ^G A2 (x 2 , Pt ) Xl G Al ( Xl ,p t ) , (1) 

is evaluated using the gluon density obtained from a simple ansatz for the unintegrated gluon distribution [9| encoding 
the saturation phenomenon: 

where x — x\ ov X2, with xi, 2 = (pt/W)e ±y ; the +(-) sign in the exponent applies to the projectile (target), and 
W = s/s is the c.m.s. energy. The factor S in cq. (TIJ is the transverse area involved in the collision (see below). The 
normalization factor K describes the conversion of partons to hadrons and is determined by a global fit to pp data at 
various energies, and to d + Au data from RHIC. 

To describe the running of QCD coupling, we use the /3-function in the one-loop approximation with Nf = 3 light 
quark flavors and Aq CD = 0.05 GeV 2 but assume that the coupling freezes at a max = 0.52 [24j ]: 



ct s (Q 2 ) = min 



12tt 



27 log^ 



(Q 2 > A| CD ). (3) 



The factor of a s (Q 2 ) in the integral ([TJ is evaluated at the scale p 2 , if this is the largest scale, or else at the lower of 
Q 2 p(y) and Q 2 S T (y). The saturation momenta are defined as 

Qliv) = Ql ^part (xo^e^y) , (4) 
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where again the_+(-) sign in the exponent applies to the projectile (target). We fix the parameters to Q = 0.6 GeV, 
xq = 0.01, and A = 0.205. In the midrapidity region of collisions at RHIC energy, this results in a gluon saturation 
momentum Q s ~ 0.68 GeV for a proton. On account of the large radius of the deuteron, we have used N part ,p = 1 
in in this case assuming that the parton substructure of the nucleon in the deuteron is not modified. For minimum 
bias d + Au collisions we multiply dN/dy by a factor of 1.52 which is our estimate for the corresponding equivalent 
number of p + Au collisions at an energy of W — 200 GeV. For pp collisions we choose the effective area S pp ~ 0.7 S p a 
somewhat smaller than for pA collisions, as suggested by the data. This may be an indication that in proton-proton 
collisions only part of the proton takes part in the interaction. On the other hand, the large nucleus makes all of the 
proton's constituents to interact. 

To evaluate the pseudo-rapidity distributions, Eq. (fTJ) needs to be rewritten using the transformation 



with the Jacobian 



, . 1, V cosh 2 ti + u 2 + sinh n 

y(v) = 7, log , ■ a — r— ( 5 ) 

z V cosn i] + u" — sinh rj 



or] Vcosh 2 77 + /i 2 



The scale p (W) is allowed to exhibit a weak energy dependence according to 



0.24 

= 0.13 + 0.32 W - 115 ' (?) 
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with W expressed in units of TcV. This parameterization reproduces approximately the "shoulder" structure of dN/ dij 
observed in symmetric pp collisions. We did not modify /x(W) for the case of pA collisions. 

The multiplicity discussed above represents the average multiplicity (iV c h) observed in collisions with a fixed number 
of participants. In experiment, the multiplicity fluctuates both due to the fluctuations in the number of participants 
and due to "intrinsic" fluctuation at fixed number of participants. To model the "intrinsic" fluctuations of the number 
of produced particles we consider the multiplicity (per unit rapidity) as a random variable distributed according to a 
negative binomial distribution, 



The quantity k which characterizes the fluctuations in the saturation approach has been estimated as be [25l [2j 

N 2 - 1 

k = K^—Q 2 s (y,W)a k (W) . (9) 

In our numerical estimates we have assumed that (Tf-(W) — <Ti n (W)/10 is proportional to the inelastic pp cross section, 
and that Q s is the saturation scale of the proton. We find that the value of k which describes best the multiplicity 
distributions in pp collisions is about k ~ 0.05. 

All observables for pA collisions finally need to be averaged also over an ensemble of N pait ^, which enters through 
eq. (HJ. We obtain the number of participants in the heavy ion target from a Monte-Carlo Glauber simulation 2 : 
assume a uniformly distributed random number < v < 1 and let 

iV partjA (6)= J2 ®(P{b-n)-u^ . (10) 

i=l—A 

Here, b is the impact parameter of the p + A collision, i.e. the transverse distance of the proton from the center of the 
target nucleus; it is a random variable with the probability density b db. The set {f^} corresponds to the coordinates of 
the nucleons in the target which are picked randomly according to a Woods-Saxon distribution. Finally, P(r) denotes 
the interaction probability of two nucleons separated by a transverse distance r; for simplicity, here we assume "hard 
sphere" nucleons: 

P(r) = e( V /^E-r) . (11) 

We use the measured values cr in (s) = 42, 52, 60, 65.75, 70.45 mb at W = 200, 900, 2360, 4400, 7000 GeV, respectively. 

Let us now present and describe our results. First we re-check the model against the RHIC data. Fig. [T] shows 
the comparison to the d Au data; in the range — 1 < r\ < 2 the agreement is satisfactory. Note that at r\ > 2 the 
saturation momentum of the projectile becomes small and so the validity of the saturation approach is questionable 
at best. Also, in the fragmentation region of the nucleus one would have to account for the contribution from valence 
quarks to improve agreement with the data. 

The centrality dependence of the charged particle multiplicity in Au+Au collisions at RHIC is shown at Fig. [H th e 
agreement is very good. The reduction of the intercept of the gluon distribution (by ~ 20% in comparison to [16J) 
allows us to reproduce well also the LHC Pb+Pb data, see Fig. O Figs. I3I4I5I show the comparison of our model 
to the pp data from the LHC on charged hadron multiplicities and multiplicity distributions at y/s = 0.9, 2.36 and 
7 TeV, respectively. The agreement is seen to be quite good. Finally, in Figs. 16171 we present our predictions for the 
upcoming p Pb run at yfs = 4.4 TcV. 

To summarize, we have presented updated predictions of the KLN model for p Pb collisions at the LHC, as well 
as comparisons to the RHIC and LHC data on hadron multiplicities and multiplicity distributions. Clearly, our 
treatment has been somewhat model-dependent and involves a few adjustable parameters. Nevertheless, our model 
does capture the emergence of a new dimensionful scale governing QCD interactions at high energies, and thus 



2 On the other hand, in AA collisions the fluctuations of TVpart do not affect the multiplicity strongly; we have calculated N pSjT t directly, 
in a "mean field approximation" , from a nuclear Woods-Saxon distribution. 
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FIG. 1: Rapidity distribution of charged particles in minimum bias d+Au collisions at W = 200 GeV. PHOBOS and BRAHMS 
data from refs. [3Q. l3lj. 




FIG. 2: Centrality dependence of the charged hadron multiplicity at rj = in AuAu collisions at W = 200 GeV [l4| and PbPb 
collisions at W = 2.76 TeV [l7|, E3] 



expresses in quantitative form the essence of the parton saturation phenomenon. The comparison of our model to the 
existing Pb Pb and the forthcoming p Pb data would also allow to deduce the amount of additional entropy produced 
during the evolution of the quark-gluon fluid in heavy ion collisions [32j . Our present treatment assumes no additional 
entropy production, which corresponds to the zero viscosity limit; a deviation from our prediction could signal the 
presence of viscous effects. 
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FIG. 3: Left: rapidity distribution of charged particles in pp collisions at W = 900 GeV. Right: Charged particle multiplicity 
distribution. ALICE and CMS data from refs. [2Sl. l29 |] . 
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FIG. 4: Left: rapidity distribution of charged particles in pp collisions at W = 2360 GeV. Right: Charged particle multiplicity 
distribution. ALICE and CMS data from refs. [H,[2S|. 
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FIG. 6: Left: rapidity distribution of charged particles in minimum bias p + Pb collisions at W = 4400 GeV. A ~ 10% overall 
normalization uncertainty is not shown explicitly. Right: Charged particle multiplicity distribution. 
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FIG. 7: Rapidity distribution of charged particles in p + Pb collisions at W = 4400 GeV for minimum bias trigger and for the 
0-50% centrality /multiplicity class. A ~ 10% overall normalization uncertainty is not shown explicitly. 
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Note added Oct. 19, 2012: 

After our paper was published, data from the ALICE Collaboration on the charged hadron multiplicity in p+Pb 
collisions at y/s = 5.02 TeV appeared [HI. The agreement between the data and our prediction is quite good over the 
entire pseudo-rapidity range of the data, —2 < r\ < +2. However, while at 77 = our prediction essentially coincides 
with the data, towards the nuclear fragmentation region, at 77 ~ —2, the prediction deviates from the data by about 
10%, see open boxes in Fig. [5] Here we would like to point out that this small discrepancy is within the uncertainty 
resulting from the Jacobian of the transformation from rapidity to pseudo-rapidity, see eqs. (|5I6I7[) . and that it can 
be easily eliminated. 
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FIG. 8: Rapidity distribution of charged particles in minimum bias p + Pb collisions at W = 5000 GeV. For these curves we 
have also accounted for the boost of the 77 = lab frame by adding a rapidity shift of Ay = | log 82/208 ~ —0.4654 to the 
right-hand-side of eq. (JSJ). mo denotes the numerator from eq. (J7J). 

In our paper, we assumed that the value of the parameter (i (determined by the typical mass and transverse 
momentum of the produced hadrons) in pA collisions is the same as in pp collisions, and is given by eq. ([7]). A more 
accurate approximation is to assume that the typical transverse mass of the produced hadrons is determined by the 
transverse momentum distribution of the produced gluons. This leads to the assumption (see e.g. eq. (27) in that 
the parameter /x 2 decreases when the saturation momentum grows, i.e. that it varies with the pseudo-rapidity and 
the participant density in the Pb nucleus as: 

O^-a^TV^-l] 
(0.13 + 0.32 IF 0115 ) 2 

where a is a parameter reflecting the rapidity dependencies of saturation momentum and of the typical mass of 
produced hadrons. The sign of the second term describes the growth of the saturation momentum of the nucleus at 
small x, towards the proton fragmentation region (positive rf). The (presumably small) magnitude of the parameter 
a is determined to large extent by non-perturbative fragmentation phenomena. In Fig. [S] we show the result of our 
computation with a = 0.0035. One can see that a small adjustment of the Jacobian of rapidity-to-pseudorapidity 
transformation that is motivated by the expected features of hadron production leads to near-perfect agreement with 
the data. 



